Numerical study of the thermodynamics of chnoatacamite 
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We study the thermodynamic properties of the clinoatacamite compound, Cu2(OH)3Cl, by consid- 
ering several approximate models. They include the Heisenberg model on (i) the uniform pyrochlore 
lattice, (ii) a very anisotropic pyrochlore lattice, and (iii) a kagome lattice weakly coupled to spins 
that sit on a triangular lattice. We utilize the exact diagonalization of small clusters with periodic 
boundary conditions and implement a numerical linked-cluster expansion approach for quantum 
lattice models with reduced symmetries, which allows us to solve model (iii) in the thermodynamic 
limit. We find a very good agreement between the experimental uniform susceptibility and the 
numerical results for models (ii) and (iii), which suggests a weak ferromagnetic coupling between 
the kagome and triangular layers in clinoatacamite. We also study thermodynamic properties in a 
geometrical transition between a planar pyrochlore lattice and the kagome lattice. 

PACS numbers: 75.10.Jm, 05.50.+q, 75.40.Cx, 05.70.-a 



I. INTRODUCTION 

The kagome and pyrochlore lattices are among the 
archetype systems for highly frustrated magnetism, with 
both lattices displaying corner-sharing frustrated plaque- 
ttes (triangles for the two-dimensional kagome lattice and 
tetrahedra for the three-dimensional pyrochlore lattice). 
There is also a geometric connection between the two 
lattices, as the pyrochlore lattice is composed of alter- 
nating kagome and triangular lattice planes stacked on 
top of each other (along the (1 1 1) body diagonal in typ- 
ical cubic spinels that display a pyrochlore lattice) . This 
leads to the possibility of structural pyrochlore lattices 
where magnetic interactions differ within kagome planes 
and between the kagome and triangular planes. 

Several three-dimensional (3D) pyrochlore lattice ma- 
terials have been shown to decouple into kagome planes 
ordered antiferromagneticallyii^ or ferromagnetically^ji 
that are fairly well isolated from the neighboring trian- 
gular plane spins. The Zn-paratacamite mineral fam- 
ily, Znj,Cu4_a:(OH)gCl2, with X > 0.3 features spin- 
i Cu^+ ions arranged on an antiferromagnetically cou- 
pled kagome lattice alternating with triangular lattice 
layers occupied by either Cu or nonmagnetic Zn ions. 
The X — 1 end member of this family, herbertsmithite, 
has attracted interest as a strong candidate to display 
a spin-liquid ground state on almost perfectly decou- 
pled two-dimensional (2D) kagome layersi^"— However, 
the best available samples are likely not stoichiometric;^ 
with a small fraction of Cu ions on the triangular lat- 
tice planes weakly (of the order of 1 K) coupled to the 
kagome planeSi^ Materials such as YBaCo407 (Ref. [lol ) 
and Yo.5Cao.5BaCo407 (Ref. [lT|) also feature alternating 
kagome and triangular layers, but with a stacking that is 
structurally distinct from the pyrochlore lattice. 

Here, we are interested in the properties of the 
mineral clinoatacamite^ a monoclinic polymorph of 
Cu2(OH)3Cl that crystallizes in the P2i/n space group 
and features spin-i Cu^"'" ions decorated on a distorted 
pyrochlore lattice. The mineral is the extension of the 



Zn-paratacamite family to a; = 0, with the monoclinic 
distortion that occurs for x < 0.3. Clinoatacamite has 
drawn attention in recent years,— ^— in part due to its 
unique pyrochlore structure and in part due to the still 
unexplained nature of successive phase transitions. Some 
studies^^'^^ have described the lattice as consisting of 
distorted kagome layers coupled weakly through trian- 
gular layers of out-of-plane spins. Others have suggested 
a pyrochlore structure with significant couplings of all 
Cu spinsji^iii Susceptibility and specific-heat measure- 
ments display two transitions upon cooling, at Tc2 = 18 K 
and Tci « 6.4 K. Long-range magnetic order-i^^i^ and a 
weak ferromagnetic moment are present below Tci- For 
temperatures Td < T < Tc2, muon oscillations are 
observedii suggesting a static local moment, which was 
originally attributed to Neel order, while neutron diffrac- 
tion experiments find no sign of ordering in this tem- 
perature range, and the specific heat anomaly at Tc2 is 
too small for the entropy change expected at an ordering 
transition. Further analysis of this unusual phase be- 
tween 6.4 and 18 K would be aided by a complete knowl- 
edge of the local bond strengths in this distorted lattice. 

In this work, we study the thermodynamic proper- 
ties of the clinoatacamite compound by considering, as 
approximate descriptions, the antiferromagnetic Heisen- 
berg model on (i) a uniform pyrochlore lattice, (ii) a 
very anisotropic pyrochlore lattice, which can be seen 
as a quasi-two-dimensional model, and (iii) a kagome 
lattice with weak ferromagnetic coupling to (otherwise 
disconnected) spins sitting on a triangular lattice, i.e., 
a two-dimensional model. We calculate the spin sus- 
ceptibility, specific heat, and entropy for these models 
using the exact diagonalization (ED) of small clusters 
with periodic boundary conditions and, only for model 
(iii), by means of an implementation of the numerical 
linked-cluster expansions (NLCEs)^'2i on an anisotropic 
checkerboard lattice that displays the required geometry. 
NLCEs yield exact results in the thermodynamic limit 
and, therefore, enable more accurate comparisons with 
experiments, while also helping us gauge finite-size ef- 



2 



fects in the exact diagonalization calculations. Using this 
method, we compare the experimental spin susceptibil- 
ity from magnetization measurements with the numerical 
results and find very good agreement in a wide range of 
temperatures. Using ED, we also examine models (i) and 
(ii) and find that results from (i) are inconsistent with ex- 
perimental data for the susceptibility. 

Furthermore, we apply the NLCE method to a 
more general anisotropic-checkerboard-lattice Heisenberg 
model, and tune the ratio of certain exchange constants 
to capture the evolution of thermodynamic quantities in 
a transition from the planar pyrochlore lattice to the 
kagome lattice. These results provide further insight on 
the nature of the spin interactions in the clinoatacamitc 
material and on the effect of frustration in the kagome 
and pyrochlore lattices. 

The paper is organized as follows: In Sec. HH we intro- 
duce the different models utilized to describe the clinoat- 
acamite compound. Section [II Al presents the pyrochlore 
lattice and its very anisotropic, quasi-two-dimensional 
version, which we use to model clinoatacamitc. Sec- 
tion IIIBI is devoted to the two-dimensional model used. 
We show how it can be seen as a Heisenberg model on 
an anisotropic checkerboard lattice, and discuss its re- 
lationship to the uniform kagome lattice and the pla- 
nar pyrochlore lattice. We also describe how NLCEs can 
be generalized to solve quantum lattice models with re- 
duced symmetries, and in particular to solve our two- 
dimensional model for clinoatacamitc. In Sec. IIIH we 
report the uniform susceptibility of clinoatacamitc as 
measured experimentally and our numerical results for 
the uniform susceptibility, specific heat, and entropy ob- 
tained within the different theoretical models by means 
of ED and/or NLCE. Finally, our results are summarized 
in Sec. |TVl 



II. APPROXIMATE MODELS FOR 
CLINOATACAMITE 

A. The isotropic and quasi-two-dimensional 
pyrochlore lattices 

Clinoatacamitc contains three crystallographically dis- 
tinct Cu sites, such that the crystal structure consists 
of kagome planes of Cu2 and Cu3 sites alternating with 
triangular planes of Cul sites>i^ These sites are distin- 
guished primarily through the Cu-O-Cu bond angle, with 
an average angle of about 96° for bonds involving a 
Cul site and an average angle of about 118° for bonds 
within the Cu2-Cu3 distorted kagome plane. (While the 
distorted lattice structure leads to some further vari- 
ation within these averages, the differences are small 
compared to the difference in average angles for the in- 
plane and between-plane cases.) On the basis of these 
differences, it has been suggested that clinoatacamitc 
should be thought of as a very- anisotropic pyrochlore 
(quasi-2D)-lattice Heisenberg model with antiferromag- 



netic kagome planes weakly coupled to triangular planes 
.JS Within this scenario, and based on bond angle consid- 
erations, the exchange interaction between layers is likely 
ferromagnetic and about one order of magnitude smaller 
than the antiferromagnctic in-plane one. 

Other works have emphasized the /i3-0H bridging ge- 
ometry of clinoatacamitc, and suggested that the mate- 
rial is best thought of as a distorted pyrochlore magnet 
with exchange interactions that are comparable in the 
kagome planes as well as between the kagome and trian- 
gular planesji^iiii In Fig. [U we show the 16-site periodic 
cluster of the pyrochlore lattice that we will use in the 
ED. 




FIG. 1. The 16-site periodic cluster of the pyrochlore lattice. 



B. The two-dimensional model 

The study of the thermodynamic properties of the 3D 
systems in Sec. Ill Al in the thermodynamic limit is very 
demanding using linked-cluster expansions. Hence, we 
will also model this material using a two-dimensional ge- 
ometry consisting of a two-layer system of kagome and 
triangular planes, as depicted in Fig. [2] For such a model, 
we can straightforwardly implement a numerical linked- 
cluster expansion, as explained below. We will show that 
this simple approximation leads to a very good agreement 
between the experimental uniform susceptibility and the 
theoretical results. 




FIG. 2. (Color online) The 16-site periodic cluster of the 
kagome lattice with extra sites inside down triangles. Pink 
(thin) bonds represent the coupling between the kagome layer 
and the sites sitting on a triangular layer in a 2D model for 
clinoatacamite. 

In order to perform a NLCE study of such a two- 
dimensional model, we start with the Heisenberg Hamil- 
tonian on the checkerboard lattice, 

i/ = ^J,,S,-S„ (1) 
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FIG. 3. (Color online) The anisotropic checkerboard lattice 
(left) and the eight realizations of the building block used in 
the square expansion NLCE (right). The shaded area rep- 
resents the kagome lattice in the limit where the red bonds 
(sides of the shaded triangles) have the same strength, J, and 
all other bonds are zero. If the strength of the blue (thick) 
bonds, J', is nonzero, and the interaction on the black (thin) 
bonds, J", is zero, then the resulting structure can represent a 
kagome lattice coupled to sites sitting on a triangular lattice. 



where Sj is the spin-i vector at site i, and is the 
strength of the exchange interaction on each bond that 
connects sites i and j. Throughout the paper, the largest 
exchange interaction in each case study sets the unit of 
energy. We consider three different types of bonds on 
the lattice, as seen in Fig. [31 There, the red (shaded) ar- 
eas make apparent the presence of an embedded kagome 
lattice in the checkerboard lattice. One can immediately 
see that by tuning the strength of the blue (thick) bonds, 
J', and black (thin) bonds, J", to zero, one captures a 
kagome lattice plus extra decoupled sites. Moreover, if 
we set J" to zero and choose J' J) to be nonzero, 
then the structure will be that of the kagome lattice cou- 
pled to sites sitting on a triangular lattice, as depicted 
in Fig. [2] Finally, if J = J' = J", then one has the pla- 
nar pyrochlore lattice. Because of the anisotropics in the 
Hamiltonian of Eq. ([T]), the usual NLCEs for the isotropic 
case^^ cannot be used here. Therefore, in the follow- 
ing, we implement a NLCE that properly deals with the 
model presented here, in which some of the symmetries 
of the lattice are broken. 



The numerical linked cluster expansion 



In linked-cluster expansionS)^^ an extensive property 
of the model per lattice site in the thermodynamic limit 
(P) is expressed in terms of contributions from all of the 
clusters, up to a certain size, that can be embedded in 
the lattice: 



recursively as the property for each cluster (V), minus 
the weights of all of its subclusters. 



P = ^L(c)7«p(c) 



(2) 



(3) 



sCc 



The contribution from each cluster (c) in Eq. ([2]) is pro- 
portional to the weight of the cluster for that property 
{wp), and to its multiplicity (L). The weight is defined 



and the multiplicity is defined as the number of ways 
that particular cluster can be embedded in the infinite 
lattice, per site. Symmetries of the lattice are often used 
in identifying topologically distinct clusters and in com- 
puting their multiplicities. This results in major sim- 
plifications of the algorithm and usually allows for ac- 
cess to larger clusters in the series. Here, we implement 
NLCEs, where V{c) is computed by means of full exact 
diagonalization,^5i3i for lattice models that break some 
of the point-group and/or translational symmetries of the 
underlying lattice. In what follows, we discuss how es- 
sentially the same expansion as for the symmetric case 
can be used for the latter cases. 

As an example, let us consider the uniform checker- 
board lattice. In the first order of the square expan- 
sion, a single crossed square has a multiplicity of 1/2 
(Refs. [2lll23h since the number of ways it can be em- 
bedded in the lattice is half the number of sites. In the 
second order, the only distinct cluster consists of two 
corner-sharing crossed squares. This cluster has a mul- 
tiplicity of 2 X 1/2, where the extra factor of two comes 
from the two possibilities for its orientation on the lat- 
tice (related by a 90° rotation), and so on.^'^ Now, con- 
sider the anisotropic lattice of Fig. [3] where, in general, 
J" ^ J' ^ J . In this case, the translational symmetries 
are reduced by a factor of two, and the point-group sym- 
metries are reduced by a factor of four, from those of the 
isotropic checkerboard lattice. So, the square expansion 
basis used for the isotropic case cannot be used for this 
lattice anymore, since the topological clusters and the 
multiplicities have changed. 

The goal is to rearrange the terms in the series to be 
able to use the square expansion basis of the isotropic 
lattice without having to redefine the topological clus- 
ters and their subclusters. Examining the problem more 
carefully reveals that the new lattice can still be tiled by 
considering two different building blocks, as opposed to 
one crossed square for the uniform lattice, which is a di- 
rect consequence of the factor-of-two reduction in trans- 
lational symmetries. These two blocks are numbered 2 
and 5 in Fig. [31 So, in the first order, one has two dis- 
tinct clusters in the sum, each with a multiplicity that is 
half of that of the single block in the first order of the 
isotropic case. This trend continues in higher orders as 
for example, in the second order, there will be four dis- 
tinct clusters, as opposed to one in the isotropic case, 
with subclusters that are the two blocks in the first or- 
der. But, just like in the first order, the multiplicities 
for each cluster are reduced by a factor proportional to 
the increasing factor in the number of clusters (four for 
the second order). Moreover, the pool of subclusters of 
these four clusters contains the same number of clusters 
of each type in the first order, namely, four from each of 
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the two building blocks. 

The above argument implies that in the expansion for 
the less symmetric checkerboard lattice, we will have dif- 
ferent realizations of clusters that existed in the expan- 
sion for the symmetric lattice, and that the latter expan- 
sion is applicable to the anisotropic case if the weight of 
each cluster is replaced by the average weight of those 
realizations. It is easy to see that the maximum num- 
ber of topologically distinct realizations of clusters in the 
isotropic square expansion for the lattice of Fig. [3] will 
be eight. This number is the same factor by which the 
point-group and translational symmetries are reduced 
from that of the isotropic checkerboard lattice. In Fig. |3J 
we have generated the eight realizations in the first order 
(among which only two are topologically distinct). Each 
of these building blocks can serve as the starting block in 
the same algorithm that generates all of the clusters in 
the expansion for the isotropic case. In fact, this guar- 
antees the generation of the eight realizations for every 
cluster in the expansion. 

The applications of this averaging scheme in NLCEs 
are not limited to the example described here. In prin- 
ciple, this method can be used in any other expansion 
(e.g., site expansion, triangle expansion, etc.), and for 
any other model with a Hamiltonian that breaks some 
symmetries of the underlying lattice. In the following 
section, we use this implementation of the NLCE method 
to calculate the properties of the lattice in Fig. [3] for val- 
ues of the exchange constant that transform its symme- 
try from a uniform planar pyrochlore lattice to near a 
kagome lattice, believed to be the appropriate model for 
the clinoatacamite compound. 



III. RESULTS 
A. Thermodynamics of clinoatacamite 

We calculate the thermodynamic properties, such as 
the specific heat, entropy, and uniform spin susceptibility 
for the Hamiltonian ([1]) on the lattice in Fig. [31 when J' = 
— O.IJ and J" = 0, to represent clinoatacamite. A 16- 
site periodic cluster of the resulting lattice is depicted in 
Fig. [2] with thick (thin) bonds representing J (J'). Note, 
however, that NLCE computes these properties directly 
for the infinite system and does not have any statistical 
or systematic errors (such as finite-size effects) within its 
region of convergence in temperature. We carry out the 
NLCE calculations to the sixth order (six building blocks 
with maximum 19 sites) of the square expansion. 

In Fig. |4l we show the spin susceptibility per site from 
the last two orders of NLCE for this system. There, we 
have also included the experimental data for this ma- 
terial. The magnetic susceptibility of a polycrystalline 
clinoatacamite sample was measured with a SQUID mag- 
netometer under an applied field of 500 Oe. The suscepti- 
bility was measured while warming from 2 to 400 K after 
field cooling. Consistent with previously published sus- 
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FIG. 4. (Color online) Uniform susceptibility per site for cli- 
noatacamite. The empty circles are the experimental results. 
The thick solid line is the last order of the NLCE for the sys- 
tem of Fig. E] with J' = -0.1 and J" = after the Wynn 
sum with two cycles of improvement.— The thick dashed line 
shows the NLCE results for the triangular expansion of the 
kagome-lattice Heisenberg model (from Ref . I2Q) ■ Thin solid 
lines are the next-to-last orders of the NLCE sums. In the 
ED for the kagome lattice, we use a 12-site cluster. For the 
case with finite J', we use the corresponding 16-site cluster 
shown in Fig. E) For the pyrochlore lattice, we use the 16-site 
cluster shown in Fig. [l] The arrow marks approximately the 
point where results from the last two orders of NLCE start 
deviating from each other. 



ceptibility results, ^'^ a weak ferromagnetic moment is ob- 
served below Tci ~ 6.4 K (not shown) and a subtle kink is 
observed in the susceptibility at Tc2 = 18 K. We will focus 
on the susceptibility above 10 K, where the experimental 
data can be compared with the numerical results. The 
experimental molar susceptibility in cgs units is related 
to the numerical one by Xexp = Cx, where the constant 
C = NAg^nl/kJ = 0.375257 J. We use J = 193 K from 
the Curie- Weiss formula, and take g = 2.14 so that the 
numerical and experimental susceptibilities match at the 
highest temperature available experimentally {T ^ 2.1). 
There is a remarkable agreement between the experiment 
and this approximate model for all of the temperatures 
above the convergence temperature of the NLCE 0.2, 
indicated by the arrow in Fig. |4]). To have a better idea 
about the effect of the extra sites of the triangular layer 
on the susceptibility of the kagome lattice, we also show 
results from a triangle-based NLCE on the kagome lattice 
with up to eight triangles i^SiSi 

It is clear that the extra sites with weak ferromag- 
netic couplings are responsible for the enhancement of 
the uniform susceptibility at low temperatures. To un- 
derstand this, we consider the limiting case where the 
sites on the triangular layer are completely decoupled 
from the ones on the kagome layer (J' = 0). In the 
thermodynamic limit, since the kagome layer contains 
only 3/4 of the sites, any property per site can be writ- 
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FIG. 5. (Color online) NLCE and ED results for (a) the spe- 
cific heat and (b) the entropy per site of the Heisenberg model 
for clinoatacamite as well as on the kagome and pyrochlore 
lattices. Lines are the same as in Fig. |4] The inset of (a) 
shows the ED results on a logarithmic temperature grid. We 
have multiplied the kagome-lattice results by | to establish a 
fair comparison with the results for the lattice of Fig. [2] and 
have included those for the 16-site pyrochlore-lattice Heisen- 
berg model, where the coupling between kagome layers is set 
to J' = —0.1. The inset in (b) is the same as the inset in (a) 
for the entropy, except that the entropy of an isolated spin is 
also properly added to that of the kagome lattice (see text). 



ten as P = jVo + f'Pkgm, where Vo is the property for 
a single site and Pkgm is the property per site for the 
kagome lattice. Therefore, in the case of susceptibility, 
a zero-temperature divergence will emerge from the sus- 
ceptibility of an isolated spin, Xo = In fact, if we 
take Xkgm to be the NLCE results for the kagome lattice 
and calculate x = -1x0 + f Xkgm, then the resulting curve 
lies very close, but slightly below, that of the NLCE with 
J' = —0.1 (see Fig. [6]), i.e., the divergence in the uni- 
form susceptibility of clinoatacamite is mostly due to the 
nearly isolated interlayer spins. However, a small nega- 
tive J' presumably produces a finite temperature order- 
ing transition in the three-dimensional material, which is 
observed in the experiments at 6 K. 

The results from ED on finite clusters with periodic 
boundary condition further support these findings. In 



Fig.m we show the spin susceptibility for the 16-site clus- 
ter of Fig. [21 and the quasi-2D model, with J' = —0.1. 
They both agree with the experimental results very well 
in the entire temperature range. We also show the ED re- 
sults for the corresponding 12-site cluster on the kagome 
lattice (which is the same cluster as in Fig. [21 but with- 
out the extra sites inside down triangles) and the uniform 
pyrochlore lattice of Fig. [1] The latter largely disagrees 
with the experimental results, invalidating the propos- 
als that clinoatacamite has such uniformity in exchange 
constants li^iii 

At this time, the lack of a nonmagnetic isostructural 
compound has made it impossible to accurately deter- 
mine the lattice contribution to the specific heat over the 
temperature range where NLCEs are valid. Therefore, we 
cannot currently compare the magnetic specific heat of 
clinoatacamite with the results of numerical calculations 
the way we have with the susceptibility. Nevertheless, 
in Fig. [SI we show the numerical results for the entropy 
and the specific heat for the models of clinoatacamite and 
the other systems discussed above, which could be use to 
compare with future experiments. Since the specific heat 
for an isolated spin is zero, the values for the J' = —0.1 
case in Fig.[5]Ja) are roughly | of those for the kagome lat- 
tice, at least for T > | J'| [see also the inset of Fig. [5ja)]. 
The position of one of the peaks, captured in the ED 
calculations both for the pure kagome and the model for 
clinoatacamite at T ^ 0.1, approximately coincides with 
the 18 K peak observed in the experiments, considering 
J ^ 193 K.— li^ii^ The existence of such a peak in the 
specific heat of the kagome-lattice Heisenberg model has 
been a topic of discussion for a long time^^i^ir— and the 
experiments with the clinoatacamite compound may have 
provided a proof of its existence. On the other hand, the 
only peak of the specific heat for the finite-size pyrochlore 
lattice from ED is at T ^ 0.3. This is inconsistent with 
the experimental results for clinoatacamite and is yet an- 
other evidence that this material is not well described by 
the uniform (or nearly uniform) pyrochlore Heisenberg 
model. 

In the inset of Fig. [SJa), we show the specific heat 
from ED on a logarithmic temperature scale and down 
to r = 0.001. The specific heat of the kagome or the py- 
rochlore lattice vanishes below T ^ 0.01, whereas a third 
peak emerges at T < 0.001 for the 2D model of clinoata- 
camite (the cluster of Fig. [2]). A similar feature also exists 
in the corresponding quasi-2D model with J' = —0.1. 
The peak moves to higher temperatures by increasing 
\J'\. Although finite-size effects often prevent ED from 
predicting, even qualitatively, the correct features of such 
models with long-range order at low temperatures, the 
appearance of this low-temperature peak due to the fi- 
nite J' may signal a possible very-low-temperature phase 
transition in the thermodynamic limit, perhaps associ- 
ated with the one observed experimentally for clinoata- 
camite at T 6 K. 

The entropies per site for the 2D and quasi-2D models 
of clinoatacamite, the kagome-lattice and the pyrochlore- 
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FIG. 6. (Color online) NLCE results for the uniform sus- 
ceptibility per site of the Heisenberg model in the transition 
between the planar pyrochlore lattice {J" = J' = J in the 
lattice of Fig. |3| and the kagome lattice with extra decoupled 
sites (J" = J' = 0). The thin dashed (dotted-dashed) line 
is the ED result for the 2D model of clinoatacamite (uniform 
pyrochlore lattice). For J' = and 0.25, black (thin) solid 
lines and color (thick) lines are the fifth and sixth orders of 
the bare sums in the expansion, respectively. For all other 
values of J', we have used Wynn extrapolation with one cy- 
cle of improvement,^ for which the thin solid and thick lines 
are the last two orders. The inset compares the uniform sus- 
ceptibility per site for the planar pyrochlore lattice and the 
kagome lattice. 



lattice Heisenberg models are shown in Fig. O^b) . Just 
like for the specific heat, we show, in the inset of Fig.[51[b), 
the low-temperature entropy of different models from the 
ED, which give us an idea of what may happen at lower 
temperatures. There, we have multiplied the entropy of 
the 12-site kagome lattice by | and added the contribu- 
tion from the isolated spins (^) to be able to properly 
compare it to the entropy of the 16-site clusters. We 
note that above T = 0.01, all entropies but the one for 
the pyrochlore lattice agree with each other. Also, as 
inferred from the specific-heat plots, a finite J' seems to 
bring about a phase transition at a very low temperature, 
after which the entropy drops to zero. 



B. Transition between planar pyrochlore and 
kagome lattices 

To gain further insights about how thermodynamic 
properties change in transitions between different frus- 
trated models, and its implications for the research on 
future materials, we study here the uniform susceptibil- 
ity, and the specific heat in the transition between the 
planar pyrochlore lattice and the kagome lattice, using 
the implementation of NLCE described in Sec. Ill Bl We 
start with the former lattice (J" — J' — J). To approach 
the kagome lattice, we simultaneously decrease J' and J" 
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FIG. 7. (Color online) NLCE results for the specific heat and 
entropy per site of the Heisenberg model on the anisotropic 
checkerboard lattice of Fig. [1 with < J" = J' < J. The 
inset compares the specific heat and the entropy per site of 
the planar pyrochlore and the kagome lattices. The lines are 
the same as in Fig. [6] 



from 1 to 0. As discussed above, the latter limit repre- 
sents the kagome-lattice Heisenberg model plus an extra 
isolated spin for every three spins in the kagome lattice, 
which is closely related to the 2D model for clinoata- 
camite. 

As can be seen in Fig. [51 the spin susceptibility of the 
planar pyrochlore lattice can even provide a good esti- 
mate for that of the 3D pyrochlore lattice (from ED), as 
the difference between the two remains relatively small 
for temperatures accessible to NLCE (T > 0.3). To show 
the proximity of the results on the other side of the tran- 
sition to the model for clinoatacamite (J' = —0.1 and 
J" — 0), we plot in Fig. [B] results from the latter from 
ED. As the spins on the triangular layer decouple from 
those on the kagome layer by decreasing J', the y di- 
vergent signature of the susceptibility of isolated spins, 
similar to what has been seen in the experiments on cli- 
noatacamite, becomes apparent. 

It is now interesting to compare the uniform suscep- 
tibility for the planar pyrochlore lattice and the pure 
kagome lattice. Within the present NLCE calculation, 
the latter can be obtained by subtracting the contribu- 
tion of isolated spins in the J" = J' = case. The results 
are shown in the inset of Fig. [51 One can clearly see there 
that the kagome lattice has a higher uniform susceptibil- 
ity than the planar pyrochlore lattice for all temperatures 
accessible within our NLCE. 

The planar pyrochlore lattice and the pure kagome lat- 
tice are two of the most frustrated lattices known. In 
Fig. [71 we show how the specific heat evolves in the tran- 
sition between them for the same parameters depicted 
in Fig. [HI Unlike the spin susceptibilities, the specific 
heat of the planar pyrochlore lattice is qualitatively dif- 
ferent from the pyrochlore lattice. In the two-dimensional 
model, as J' and J" decrease, the high-temperature peak 
is suppressed. However, this is largely due to the fact 
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that one-fourth of the spins in the system are decoupled 
from the lattice in the limit of J" = J' = and, there- 
fore, have vanishing specific heat. Consequently, if one 
compares the entropy and specific heat per site of the 
planar-pyrochlore-lattice and the kagome-lattice Heisen- 
berg models (inset in Fig. [7]) , one sees that their values 
are in fact very close for all of the temperatures calcu- 
lated here. Interestingly, this shows that both lattices 
have a very similar degree of frustration. 

IV. CONCLUSIONS 

We have presented a numerical study of the thermody- 
namic properties for models of the clinoatacamite com- 
pound. In particular, we computed the spin suscepti- 
bility, entropy, and specific heat, using the ED of finite 
periodic clusters and an implementation of the NLCEs 
that properly deals with the breaking of lattice symme- 
tries introduced by the particular model Hamiltonian of 
interest. We find an excellent agreement between the ex- 
perimental uniform susceptibility of clinoatacamite from 



magnetic measurements and our numerical results for the 
Heisenberg model on a lattice that consists of a kagome 
layer, coupled weakly to a triangular layer. Together with 
a study of the entropy and the specific heat of the kagome 
and pyrochlore lattices, we provide strong evidence that 
clinoatacamite has a pyrochlore structure with only weak 
ferromagnetic coupling between its kagome layers. Em- 
ploying our generalized NLCE, we also studied the above 
thermodynamic quantities in a transition between the 
planar pyrochlore lattice, which has a uniform suscepti- 
bility similar to that of the pyrochlore lattice, and the 
kagome lattice plus isolated spins, closely related to the 
model for the clinoatacamite compound. 
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